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Abstract 

We continue the study of a model for heat conduction |^ consist- 
ing of a chain of non-hnear osciUators coupled to two Hamiltonian 
heat reservoirs at different temperatures. We establish existence of a 
Liapunov function for the chain dynamics and use it to show exponen- 
tially fast convergence of the dynamics to a unique stationary state. 
Ingredients of the proof are the reduction of the infinite dimensional 
dynamics to a finite-dimensional stochastic process as well as a bound 
on the propagation of energy in chains of anharmonic oscillators. 

1 Introduction 

In its present state, non-equilibrium statistical mechanics is lacking the firm 
theoretical foundations that equilibrium statistical mechanics has. This is 
due, perhaps, to the extremely great variety of physical phenomena that 
non-equilibrium statistical mechanics describes. We will concentrate here 
on a system which is maintained, by suitable forces, in a state far from 
equilibrium. In such an idealization, the non-equilibrium phenomena, can be 
described by stationary non-equilibrium states (SNS), which are the analog 
of canonical or microcanonical states of equilibrium. 

Recently many works have been devoted to the rigorous study of SNS. 
Two main streams are emerging. In the first approach, for open systems^ a 
system is driven out of equilibrium by interacting with several reservoirs at 
different temperatures. In the second approach, for thermostated systems , 
a system is driven out of equilibrium by non-Hamiltonian forces and con- 
strained to a compact energy surface by Gaussian (or others) thermostats 
P, p^]. One should view both approaches as two different idealizations of 
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the same physical situation, in the same spirit as the equivalence of ensem- 
bles in equilibrium statistical mechanics. But for the moment, the extent to 
which both approaches are equivalent remains a largely open problem. 

We consider here an open system, a model of heat conduction consist- 
ing of a finite-dimensional classical Hamiltonian model, a one-dimensional 
finite lattice of anharmonic oscillators (referred to as the chain), coupled, at 
the boundaries only, to two reservoirs of classical non-interacting phonons 
at positive and different temperatures. We believe this model to be quite 
realistic, in particular it is completely Hamiltonian and non-linear. 



2§). First 
Sim- 



This model goes back (in the linear case) to |^] (see also [22 
rigorous results for anharmonic models appear Q and further in 
ilar models in classical and quantum mechanics have attracted attention in 
the last few years, mostly for systems coupled to a single reservoir at zero 
or positive temperature, i.e., for systems near thermal equilibrium (see e.g. 



[Il2| , |13| , y, [15|, In our case, with two reservoirs, no Gibbs Ansatz is avail- 
able and in general, even the very existence of a (non-equilibrium) stationary 
state is a mathematically challenging question which requires a sufficiently 
deep understanding of the dynamics. For the model at hand, conditions for 
the existence of the SNS have been given in Q and generalized in The 
uniqueness of the SNS as well as the strict positivity of entropy production 
(or heat flux) have been proved in Q. The leading asymptotics of the in- 



variant measure (for low temperatures) are studied in |21] and shown to be 
described by a variational principle. 

Under suitable assumptions on the chain interactions and its interactions 
with the reservoirs, we establish the existence of a Liapunov function for the 
chain dynamics. We then use this Liapunov function to establish that the 
relaxation to the SNS occurs at an exponential rate, and finally we prove 
that the system has a spectral gap (using probabilistic techniques developed 



by Meyn and Tweedie in [18|). 



The Hamiltonian of the model has the form 

H = Hb + Hs + Hj. (1) 

The two reservoirs of free phonons are described by wave equations in R'^ 
with Hamiltonian 

Hb = H{ipl,ttl) + H{ipr,ttr) , 
Hiv,7T) = iy'^ix(|V(^(x)|2 + |vr(x)|2), 

where L and R stand for the "left" and "right" reservoirs, respectively. The 
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Hamiltonian describing the chain of length n is given by 

n 2 

1=1 ^ 

n n—1 
j=l i=l 

where (pi, gj) G R'^ x R'^ are the coordinates and momenta of the i^^ particle 
of the chain. The phase space of the chain is R^'^". The interaction between 
the chain and the reservoirs occurs at the boundaries only and is of dipole- 
type 

Hi = qi- J dxV(pL{x)pL{x) + qn- j dx\/ipR{x)pR{x) 

where pi and pR are coupling functions ("charge densities") which we will 
assume spherically symmetric. 

Our assumptions on the anharmonic lattice described by Hs{p,q) are 
the following: 

• HI Growth at infinity: The potentials U^^\x) and U^'^\x) are 
and grow at infinity like and ||x||^2: There exist constants Cj, 

Di, i = 1,2 such that 

lim A-^»C/(*)(Ax) = a»||xf\ (2) 

A— >oo 

lim A-'=»+^V?7W(Ax) = a^'hi\\x\\^^-^x , (3) 



A— >oo 



||52[7»(x)|| < {Q + DiV{x))' (4) 

where || • || in Eq. (^) denotes some matrix- norm. 
Moreover we will assume that 

k2 > ki > 2, 

so that, for large ||x|| the interaction potential C/^^^ is "stiffer" than 
the one-body potential U^^\ It follows from Eqs. (||) and (^) that the 
critical set of V{q), i.e, the set {q : W{q) = 0} is a compact set. 

H2 Non-degeneracy: The coupling potential between nearest neigh- 
bors is non-degenerate in the following sense. For x G R'^ and 
m = 1,2, • • •, let A^"^\x) : R'^ R'^'" denote the linear maps given 
by 

d Qm+1^{2) 
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We assume that for each x G R"^ there exists rriQ such that 
Rank(^(i)(x),---^(™")(x)) = d. 



The class of couphng functions pi, i £ {L,R} we can ahow is relatively 
restrictive: 

• H3 Rationality of the coupling: Let pi denote the Fourier trans- 
form of Pi. We assume that 



where Qi, i € {L,R} are polynomials with real coefficients and no 
roots on the real axis. In particular, if /cq is a root of Qi, then so are 
—ko, /cq and — /cq. 

Under these conditions we have the following result (a more detailed 
and precise statement will be given in the next section). Let F{p,q) be 
an observable on the phase space of the chain, for example any function 
with at most polynomial growth (no smoothness is required). We denote as 
{p{t), q{t)) the solution of the Hamiltonian equation of motion with Hamil- 
tonian (|l|) and initial conditions {p, q). Of course (p(t), q{t)) depends also on 
the variable of the reservoirs, though only through their initial conditions 
{t^Li^Li'^Ri'^r)- We introduce the temperature by making the assumption 
that the initial conditions of the reservoirs are distributed according to ther- 
mal equilibrium at temperature Tr and Tl respectively and we denote {■)lr 
as the corresponding average. 

Theorem 1.1 Under Conditions HI— H3, there is a measure ^'{dp, dq) with 
a smooth everywhere positive density such that the Law of Large Numbers 
holds: 



for almost all initial conditions {ttl, tt/j, (Pr) of the reservoirs and for all 
initial conditions {p, q) of the chain. Moreover there exist a constant r > 1 
and a function C (p, q) with f C dv < oo such that 



for all initial conditions {p, q) . That is, if we average over the initial condi- 
tions of the reservoirs the convergence is exponential. 



1 
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Note that the ergodic properties stated in Theorem |1.1| hold not only for 
i/-almost every initial condition {p,q), but in fact for every {p,q). 

The existence of a (unique) stationary state was proved for (exactly 
solvable) quadratic harmonic potentials V{q) in [^], for fci = /c2 = 2 (i.e., 
for potential which are quadratic at infinity) in |^ and generalized to the 
case /c2 > fci > 2 in What is really new here is that we prove that 
the convergence occurs exponentially fast and we also weaken slightly the 
conditions on the potential (in particular the case ki = k2 is allowed and 
our condition H2 on C/^^^ is weaker than the one used in |, 0, |]). Our 
methods also differ notably from those used in ^ ; in fact we reprove the 
existence of the SNS (with a shorter and more constructive proof than in 

^]) and, at the same time, we prove much stronger ergodic properties. 

We devote the rest of this section to a brief discussion of the assumptions 
HI — H3. Since the reservoirs are free phonons gases and since we make a 
statistical assumption on the initial condition of the reservoirs, one can inte- 
grate out the variables of the reservoirs yielding random integro-differential 
equations for the variable {p, q). Our assumption H3 of rational coupling is, 
in effect, a Markovian assumption: with such coupling one can eliminate the 
memory terms by adding a finite number of auxiliary variables to obtain a 
system of Markovian stochastic differential equations on the extended phase 
space consisting of the dynamical variables {p, q) together with the auxiliary 
variables. The main (new) ingredient in our proof is then the construction 
of a Liapunov function for the system, which implies, using probabilistic 
methods developed in |2^, 18|, the exponential convergence. 

To explain the construction of a Liapunov function, note that the dynam- 
ics of the chain in the bulk is simply Hamiltonian, while at the boundaries 
the action of the reservoirs results into two distinct forces. There are dis- 
sipative forces which correspond to the fact that the energy of the chain 
dissipates into the reservoirs. This force is independent of the temperature. 
On the other hand since the reservoirs are infinite and at positive tempera- 
tures, they exert (random) forces at the boundaries of the chain and these 
forces turn out to be proportional to the temperatures of the reservoirs. 

The construction of the Liapunov function proceeds in two steps. In a 
first step we neglect completely the random force, only dissipation acts. This 
corresponds to dynamics at temperature zero, and one can prove that the 
energy decreases and that the system relaxes to a (local) equilibrium of the 
Hamiltonian H{p,q). We establish the rate at which this relaxation takes 
place (at sufficiently high energies). In the second step we consider the com- 
plete dynamics and we show that for energies which are much higher than 
the temperatures of the reservoirs, the random force is essentially negligible 
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with respect to the dissipation. This means that except for (exponentially) 
rare excursions the system spends most of its time in a compact neighbor- 
hood of the equilibrium points. On the other hand, in this compact set, i.e., 
at energies of order of the temperatures of the reservoirs, the dynamics is 
essentially determined by the fluctuations and to prove exponential conver- 
gence to a SNS one has to show that the fluctuations are such that every 
part of the phase space is visited by the dynamics. To summarize, we control 
the dynamics at any temperature by the dynamics at zero temperature. 

This allows one to understand the meaning of our assumptions on the 
potential V{q). If we suppose that the energy has an infinite number of local 
minima tending to infinity, the zero temperature (long-time) dynamics is not 
confined to a compact energy domain and our argument fails. With regard 
to the condition k2 > ki in HI on the exponents of the potentials, since 
the results of and the rigorous proofs of ^, it is known that stable 
(in the sense of Nekhoroshev) localized states exist in non-linear lattices. 
Consider, for example, an infinite chain of oscillators (without reservoirs). 
Numerically and in certain cases rigorously [17|, one can show the existence 
of breathers, i.e., of solutions which are spatially (exponentially) localized 
and time-periodic. Although the breathers occurs both for ki > k2 and 
k2 > ki they behave differently at high energies. For ki > k2, the higher 
the energy, the more localized the breather get (hard breathers), while for 
k2 > ki, as the the energy gets bigger the breather become less and less 
localized (soft breathers). In fact a key point of our analysis is to show 
that at high energy, if the energy E of the initial condition is localized away 
from the boundary, then after a time of order one, the oscillators at the 
boundaries carry at least an energy of order E"^^^^ so that the chain system 
energy can relax into the reservoirs. 

Although we believe that the existence of a SNS probably may not de- 
pend too much on these localization phenomena, the rate of convergence to 
the SNS presumably does. Our approach of controlling the dynamics by the 
the zero-temperature dynamics may not be adequate if the condition HI 
fails to hold and so more refined estimates on the dynamics are needed to 
show that these localized states might be in fact destroyed by the coupling 
to the reservoirs. 

As regards the organization of this paper, Sec. |^ presents the effective 
stochastic differential equations for the chain, a discussion of allowable in- 
teractions between the reservoirs and the chain and a concise statement. 
Theorem 2A, of the exponential convergence. In Sec. |3| we discuss the dis- 
sipative deterministic system (corresponding to reservoirs at temperature 
0) , Theorem |3.3| , and then we show the extent to which the random paths 
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follow the deterministic ones, Propositio n |3.7| . We give a lower bound on 
the random energy dissipation, Corollary |3.8| . We then conclude Sec. ^ by 
providing the Liapunov function, Theorem 3.10| , and bounds on the expo- 
nential hitting times on (sufficiently large) compact sets, Theorem 0.11 . In 
Sec. ^ we prove the random process has a smooth law and at most one 
ergodic component, improving slightly results of 0, ^. Finally in Sec. 
^ we conclude the proof of Theorem 2T by invoking results of on the 
ergodic theory of Markov process. 



2 Effective Equations 

We first give a precise description of the reservoirs and of their coupling 
to the system and derive the stochastic equations which we will study. A 
free phonon gas is described by a linear wave equation in R'^, i.e., by the 
pair of real fields (j){x) = {ip{x),7r{x)), x G R'^. We define the norm 

ll'/'lP = J dx \\/(j){x)\'^ + |'/r(x)p and denote (•, •) the corresponding scalar 
product. The phase space of the reservoirs at finite energy is the real Hilbert 
space of functions (j){x) such that the energy Hb{4)) = ||^|p/2 is finite and 
the equations of motion are 

<p{t,x) = C(l){t,x), C = 

In order to describe the coupling of the reservoir to the system, let us con- 
sider first a single confined particle in R'^ with Hamiltonian Hs{p,q) = 
p'^/2-\-V{q). As the Hamiltonian for the coupled system particle plus reser- 
voirs, we have 

H{4>,p,q) = ^Uf + V{q)+q- J dxVip{x)p{x) 

= HB{cl))+Hsip,q)+q-{cP,a), 

where p{x) is a real rotation invariant function and a = {a^^\ • • • , a^^^) is, 
in Fourier space, given by 

We introduce the covariance matrix = (exp (£t)a(*) , a^^)). A simple 

computation shows that 
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and we define a coupling constant A by putting = C^^^\0) = J dk\p{k)\'^. 
The equations of motion of the coupled system are 



m = p{t), 

m = -vy(g(t))-(0,a), 
4>{t,k) = C{(b{t,k)+q{t)-a{k)) . (5) 

With the change of variables ilj{k) = 4>{k) + q ■ a{k), Eqs. (|5|) become 

m = p{t), 

p(t) = -Vyeff(?(i))-(V',a), 
ii:{t,k) = Cip{t,k)+p{t)-a{k), (6) 

where Vcsiq) = ViQ) — A^g^/2. Integrating the last of Eqs. (^) with initial 
condition Tpo{k) one finds 

il^{t,k) = e^Vo(^)+ t dse'^^*-'^a{k)-p{s). 
Jo 

and inserting into the second of Eqs. (^) gives 

m = Pit), 

p[t) = -VV,siq{t)) - [ dsC{t - s)p{s) - {iPo, e-'^'a) . (7) 



If we now assume that, at time t = 0, the reservoir is at temperature T, 
then ipQ is distributed according to the Gaussian measure with covariance 
T(-,-) and then ^{t) = {tpo, e~'-'^a) is a d-dimensional stationary Gaussian 
process with mean and covariance C{t — s). Note that the covariance 
itself appears in the deterministic memory term on the r.h.s. of Eq.(^) 
(fluctuation-dissipation relation) . 

By assumption H3 there is a polynomial p{u) which is a real function of 
iu and which has its roots in the lower half plane such that 

/OO 1 

Note that this is a Markovian assumption Q: ^(t) is Markovian in the sense 
that we have the identity p[—id/ dt)^[t) = uj{t), where co{t) is a white noise, 
i.e, the joint motion of d^^(t) / df^ , < m < degp — 1 is a (Gaussian) 
Markov process. This assumption together with the fluctuation-dissipation 
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relation permits, by extending the phase space with a finite number of vari- 
ables, to rewrite the integro-difFerential equations (0) as a Markov process. 
Note that ^(t) can be written as 

/oo r 
k{t-t')(Lj{t'), k{t) = due'''^p{u)-^ . 
-oo J 

with k{t) = for t < 0. For example \i p{u) oc iu + j then C{t) = A^e"'*''*'. 
Introducing the variable r defined by 

Xr{t) = f ds C{t - s)p{s) + /* k{t - t')doj{t') , 

JO J-oo 

then we obtain from Eqs.(|^ the set of Markovian differential equations: 

m = Pit), 

p{t) = -VV,s{q{t)) - Xr{t) , 

dr{t) = {-jr{t) + Xp{t))dt + {2TjY/'^duj{t). (8) 

If p{u) oc {iu + 7 + ia){iu + 7 — icr) then C{t) = \^ cos{at)e~"'^*^ and intro- 
ducing the two auxiliary variables r and s defined by 

Xr{t) = A2 /* ds cos{a{t - s))e-^^^-'^p{s) 
Jo 

+ (rA27)i/2 /■* cos(cT(t - s))e~Tl*-^l(iw(s) , 
J —00 

Xs{t) = A2 /* dt sm{a{t - s))e-'^\^-'\p{s) 
Jo 

+ (rA27)i/2 /■* dt sm{a{t - s))e-^\'~'Uuis) , 
J —00 

we obtain then the set of Markovian differential equations: 

m = p{t), 

Pit) = -VV,s{q{t)) - Xrit) , 

dr{t) = {--ir{t)-as{t) + Xp{t))dt + {2T-f)^/'^duo{t), 

s{t) = -js{t) + ar{t) . (9) 

Obviously others similar set of equations can be derived for arbitrary poly- 
nomial p{u). 

Another coupling which we could easily handle with our methods oc- 
curs in the following limiting case, see |^]. Formally one wants to take 
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C{t) = r]^6{t). Note that this corresponds to a couphng function with 
|;9(fc)P = 1 in which case = cxd. A possible hmiting procedure consists 
in taking a sequence of covariances tending to a delta function and at the 
same time suitably rescaling the coupling (see Q). In this case one obtains 
the Langevin equations which serve as commonly-used model system with 
reservoir in the physics literature, 

m = p{t), 

dp{t) = {-VVMt))-v'p{t))dt + {2Tr,^fl^du:{t). (10) 

The derivation of the effective equations for the chain is a straightforward 
generalization of the above computations. Our techniques apply equally well 
to any of the couplings above. However, for simplicity, we will only consider 
the case where the couplings to both reservoirs satisfy \pi{k)\^ oc /c^ + 7^, 
i = L,R. For notational simplicity we set Ti = Tl and r„ = Tr, we 
denote ri and r„ as the two auxiliary variables and we will use the notations 
r = (ri,r„), and x = {p,q,r) G X = R2c({n+i)_ jj^ ^j^jg ^g^gg obtain the 
set of Markovian stochastic differential equations given by 

^1 = Pi, 

Pi = -"^qiVcsiq) - >^ri , 

dn = (-7ri + Xpi) dt + {2Ti-if/^duJi , 

qj = Pj, j = 2, . . . , n - 1 , 

Pj = -Vq,Kfr(g) i = 2, . . . , n - 1 , 

— Pn 1 

Pn = -Vg„V;ff(q') - Ar„ , 

drn = {-^rn + \pn)dt + {2Tn^f'^dujn, (11) 

where Vesiq) = V{q) - \^ql/2 - \^ql/2. From now on, for notational sim- 
plicity we will suppress the index "eff " and consider V = T4ff as our potential 
energy. 

It will be useful to introduce the following notation. We define the linear 
maps A : R'^" B?'^ by A(2;i, . . . , x^) = (Axi, Ax„) and T : R^'^ R^'^ by 



T{x,y) = (Tix,T„y). With this we can rewrite Eq.(ll) in the compact form 

q = p, 

p = -VqV - A^r , 

dr = {--ir + Kp)dt + {2'yTf/'^doj. (12) 
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The solution x{t) of Eq.(|T|) is a Markov process. We denote T* as the 
associated semigroup, 

T'fix) = E,[/(x(t)], 

with generator 

L = 7 {VrTVr - rVr) + {ApVr - rAVp) + {pVg - {VqV{q))Vp) , (13) 

and Pt{x, dy) as the transition probabihty of the Markov process x{t). There 
is a natural energy function which is associated to Eq.(^), given by 

Gip,q,r) = y + H{p,q) . 
A straightforward computation shows that in the special case Ti = T„ = T 

^-lg-G(p,g,r)/T 

is an invariant measure for the Markov process x{t). 

Given a function W : X ^ TL satisfying 1^ > 1 we consider the following 
weighted total variation norm || • \\w given by 



|7r||iy = sup 
\f\<w 



J fdn\ , (14) 



for any (signed) measure ir. We introduce norms || • We and Banach spaces 
Lf{X) given by 

11/11, = sup L^(X) = {/: 11/11, <oo}, (15) 

and write \\K\\0 for the norm of an operator K : L'^{X) L^{X). 



Theorem 1.1 is a direct consequence of the following result 



Theorem 2.1 Assume that HI and H2 hold. The Markov process x{t) 
which solves has smooth transition probability densities, Pt{x,dy) = 

pt{x,y)dy, with pt{x,y) £ C°°{{0,oo)x X x X). The Markov process x{t) has 
a unique invariant measure /U, and fi has a C°° everywhere positive density. 
For any 9 with < 9 < (max{Ti, T„})~^ there exist constants r = r{9) > 1 
and R = R{9) < oo such that 

\\Pt{x,-)-f^\\e.p(9G) < Rr~'e^p{9G{x)), (16) 
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for all X G X, (exponential convergence to the SNS) or equivalently 

(spectral gap). Furthermore for all functions f , g with g'^ G L'^{X) and 
all t > we have 



gT^f dn- J fdn J gdfj. 



tu .2||1/2|, 2||l/2 



< Rr-^WfYe'llg 



(exponential decay of correlations in the SNS). 

The convergence in the weighted variation norm, Eq. (|l^), imphes that 



the Law of large Numbers holds |1C, 18 1. 



Corollary 2.2 Under assumptions HI and H2 x{t) satisfies the Law of 
Large Numbers: For all initial conditions x G X and all f G L^ {X, dfx) 



hm 1 r f{x{t))dt = ffdi^L 
T^oo 1 Jo J 



almost surely. 



The convergence of the transition probabilities as given in (^) is shown 
in 1 15] to follow from the following properties: 



Strong Feller property The diffusion process is strong Feller, i.e., 
the semigroup T* maps bounded measurable functions into continuous 
functions. 



This is a consequence of the hypoellipticity of the diffusion x{t), which fol- 
lows from Condition H2, see Section ^. 

• Small-time open set accessibility. For all t > 0, all x G X and all 
open set A C X we have Pt{x, A) > 0. 

This means that the Markov process is "strongly aperiodic" . In particular, 
combined with the strong Feller property it implies uniqueness of the in- 
variant measure. This property is discussed in Section Q using the support 
theorem of and explicit computations. This generalizes (slightly) the 
result obtained in 0]. 
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• Liapunov function and hitting times Fix s > arbitrary. Set 
W = exp(6'G) and choose 6 with <e < (max{ri, r„})~^ Then W 
is a Liapunov function for the Markov chain {x{ns)}n>o- W > 1, W 
has compact level sets and there is a compact set U, (depending on s 
and 9) and constants k < 1 and b < oo, (both depending on U, s and 
9) such that 

T'W{x) < KW{x) + blu{x), (17) 
where Ijj denotes the indicator function of the set U. 

The existence of a Liapunov function is the main technical result of this 
paper (see Section ^) and the condition HI is crucial to obtain it. Note that 
the time derivative of the (averaged) energy 

jE^AGixm = 7E.[TV(T)-r2(t)], 

is not necessarily negative. But it is the follows from our analysis 

below that, that for t > 0, E^[G(x(t))-G(x)] < -cG{xf/^^ for x sufficiently 
large. 

A nice interpretation of a Liapunov bound of the form ( p^ ) is in terms of 
hitting times. Let tjj denote the first time the diffusion x{t) hits the set U; 
then Eq. ( p^ implies that tjj is exponentially bounded. We will show that 
for any a > 0, no matter how large, we can find a compact set U = U{a) 
such that 

E,.[e'^-^] < cxD, 

for all X £ X. So except for exponentially rare excursions the Markov pro- 
cess x{t) lives on the compact set U. Combined with the fact that the pro- 
cess has a smooth law, this provides an intuitive picture of the exponential 



convergence result of Theorem 2.1 



3 Liapunov Function and Hitting Times 

3.1 Scaling and Deterministic Energy Dissipation 

We first consider the question of energy dissipation for the following deter- 
ministic equations 

q = P, 

p = -V,V{q)-A^r, 

r = —^r + Ap, (18) 
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obtained from Eq.(0) by setting Ti = T„ = 0, corresponding an initial 
condition of the reservoirs with energy 0. A simple computation shows that 
the energy G{p,q,r) is non-increasing along the flow x{t) = (p{t), q{t),r{t)) 
given by Eq.(|l8|): 

j^G{pit),qit),r{t)) = -^r\t) < 0. 

We now show by a scaling argument that for any initial condition with 
sufficiently high energy, after a small time, a substantial amount of energy 
is dissipated. 

At high energy, the two-body interaction U^'^^ in the potential domi- 
nates the term C/^^^ since /c2 ^ and so for an initial condition with en- 
ergy G{x) = E, the natural time scale - essentially the period of a single 
one-dimensional oscillator in the potential \q\^^ - is _E^/'^2~i/2. "We scale a 
solution of Eq.(p!8[) with initial energy E as follows 

p{t) = E-^p{E^~h) , 
q{t) = E'^q{E'^~h) , 

= E~^r{E^'h) . (19) 
Accordingly the energy scales as G{p,q,r) = EGe{p, q,f), where 

GEip,q,r) = Ei-'^- + ^ + VEiq), 

n n—1 
i=l 1=1 

U(^{x) = E-^U^'^E^x), i = l,2. 
The equations of motion for the rescaled variables are 
q = P, 

P = -VgVE{q)-E^2'^A^r, 
1 

r = -£''=2 27f + Ap. (20) 
By assumption HI, as ^ oo the rescaled energy becomes 
Goo{p,q,r) = lim GE{p,q,f) 

E^oo 




fjl + Voo{q) fci = ^2 > 2 or fc2 > fci > 2 

fV2 + pV2 + V;o(g) ki = k2 = 2 
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where 



Ea^^^\\q^-qi+l\\^' k2>h>2 



The equations of motion scale in this hmit to 



q = P, 

p = -VgFoo®, 

f = Ap, (21) 



in the case k2 > 2, while they scale to 

q = P, 



in the case ki = k2 = 2. 



P = -VgVoo{q)-A^r, 

f = -jr + Ap , (22) 



Remark 3.1 The scaling for the p and q is natural due to the Hamiltonian 
nature of the problem, but the scaling of r has a certain amount of arbi- 
trariness. Since G is quadratic in r, it might appear natural to scale r with 
a factor E~^^'^ instead of E~^^^'^ as we do. On the other hand, the very 
definition of r as an integral of p suggests that r should scale as q, as we 
have chosen. 



Remark 3.2 Had we supposed, instead of HI, that ki > k2, then the natu- 
ral time scale at high energy would be E^^/'^i^^/^. Scaling the variables (with 
k2 replaced by ki would yield the limiting Hamiltonian p'^/2 + J2 (^^^M\qi\\^^ i 
i.e., the Hamiltonian of n uncoupled oscillators. So in this case, at high en- 
ergy, essentially no energy is transmitted through the chain. While this does 
not necessary preclude the existence of an invariant measure, we expect in 
this case the convergence to a SNS to be much slower. In any case even the 
existence of the SNS in this case remains an open problem. 



Theorem 3.3 Given r > fixed there are constants c > and Eq < 00 
such that for any x with G{x) = E > Eq and any solution x{t) of Eg. ^TB^ 
with x{0) = X we have the estimate, for tE = E^^^^'^^'^t, 

G{x{tE)) -E < -cE^~^ . (23) 
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Remark 3.4 In view of Eq. (p3D, this shows that r is at least typicahy 
O(^Vfe) on the time interval [0, ^i/*=2-i/2^]^ 

Proof: Given a solution of Eq.(|l^) with initial condition x of energy G{x) = 
E, we use the scaling given by Eq.(Q9h and we obtain 



G{x{tE))-E = -7 dtr'^it) = -7^4 2 f dtf'^{t), (24) 

JO JO 



where f{t) is the solution of Eq.(pOD with initial condition x of (rescaled) 
energy Ge{x) = 1. By Assumption H2 we may choose Eq so large that for 
E > Eq the critical points of Ge are contained in, say, the set {Ge < 1/2}- 
For a fixed E and x with G{x) = E, we show that there is a constant 
Cx,E > such that 

'^tf^t) > c^^E- (25) 




The proof is by contradiction, c.f. |21|. Suppose that Jq dt f (t) = 0, then we 



have r{t) = 0, for all t G [0, r]. Prom the third equation in ( ^0[ ) we conclude 
that pi{t) = Pn{t) = for all t G [0, r], and so from the first equation in ( |20| ) 
we see that qi{t) and qn{t) are constant on [0,r]. The second equation in 
(pOD gives then 



= f>^{t) = -V^.Vim) = -V,-i?7«(gi(t))-V,-,f/(2)(^~i(t)-g2(t)), 

together with a similar equation for pn- By our assumption HI the map 
V[/(2) has a right inverse g locally bounded and measurable and thus we 
obtain 

ut) = m - giu^'Hrn)) ■ 

Since qi is constant, this implies that ^2 is also constant on [0, r]. Similarly 
we see that qn-i is constant on [0, r]. Using again the first equation in 
(pO|) we obtain now P2{t) = Pn-i{t) = for all t G [0,r]. Inductively one 
concludes that f = implies p = and VqV = and thus the initial 
condition x is a critical point of G^;. This contradicts our assumption and 
Eq. dl) follows. 

Now for given E, the energy surface G_e is compact. Using the continuity 
of the solutions of O.D.E with respect to initial conditions we conclude that 
there is a constant ce > such that 



inf / dtf'^it) > CE 
Gk=i}Jq 



xe{Gi3=i} 
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Finally we investigate the dependence on E ce- We note that for E = oo, 



Goo lias a well-defined limit given by Eq.(21) and the rescaled equations of 
motion, in the limit — > oo, are given by Eqs. (^Tj) in the case /c2 > 2 and 
by Eq. (|2^ ) in the case ki = k2 = 2. Except in the case ki = k2 = 2 the 
energy surface {Goo = 1} is not compact. However, in the case ki = k2> 2, 
the Hamiltonian Goo and the equation of motion are invariant under the 
translation r ^ r + a, for any a G R^'^. And in the case k2 > ki > 2 
the Hamiltonian Goo and the equation of motion are invariant under the 
translation ri-^r + aq^q + b, for any a G R^*^ and b £ R'^". The quotient 
of the energy surface {Goo = 1} by these translation, is compact. 

Note that for a given x £ {Goo = 1} a similar argument as above show 
that dt{r + a)"^ > 0, for any a > and since this integral clearly goes to 
oo as a — > oo there exists a constant Coo > such that 



inf [ f'^{t) dt > Co 
Goo=iWO 



ie{Goo=l} 



Using again that the solution of O.D.E depends smoothly on its parameters, 
we obtain ^ 

inf inf / dtf^{t) > c. 



This estimate, together with Eq. (|24| ) gives the conclusion of Theorem 3^ . 



3.2 Approximate Deterministic Behavior of Random Paths 

In this section we show, that at sufficiently high energies, the overwhelming 
majority of the random paths x{t) = x{t,u>) solving Eqs.(|l2|) follows very 
closely the deterministic paths Xdet solving Eqs.(|l^). As a consequence, 
for most random paths the same amount of energy is dissipated into the 
reservoirs as for the corresponding deterministic ones. We need the following 
a priori "no-runaway" bound on the growth of G{x{t)). 

Lemma 3.5 Let 6 < (max{ri, Then E^[exp ((9G(x(t)))] is well- 

defined and satisfies the bound 

Bx[exp{eG{x{t)))] < exp(7Tr(r)et)exp(0G(x)). (26) 

Moreover for any x with G{x) = E and any 6 > we have the estimate 

P^. J sup G{x{s)) >{1 + 5)e\ < exp (7TV(r)0t) exp (-SOE) . (27) 

lo<s<t I 
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Remark 3.6 The lemma shows that for E sufficiently large, with very high 
probability, G{x{t)) = 0{E) if G{x) = E. The assumption on 9 here arises 



naturally in the proof, where we need (1 — 6T) > 0, cf. Eq. (25). 



Proof: For 6 < (max{Ti, T„}) ^ we have the bound (the generator L is given 
by Eq. dlD) 

Lexp{eG{x)) = -fO e^p {eG{x)) {Tv{T) - r{l - eT)r) 

< j6Tv{T) exp {9G{x)) , (28) 

so that for the function W{t, x) = exp {—^9Tr{T)t) exp {9G{x)) we have the 
inequality [dt + L)W{t,x) < 0. We denote aR as the exit time from the 
set {G{x) < R}, i.e., an = in{{t > 0,G{x{t)) > R}. If the initial condition 
X satisfies G{x) = E < R, we denote XR(t) the process which is stopped 
when it exits {G{x) < R}, i.e., XR^t) = x{t) for t < aji and XR{t) = x{aii) 
for t > CT/j. We set (TR{t) = mm{aji,t} and applying Ito's formula with 
stopping time to the function W{t, x) we obtain 

[exp {9Gix{aR{tm exp {-j9TT{T)aR{t))] - exp {9G{x)) < , 

thus 

[exp {9G{x{aR{t))))] < exp (7^Tr(r)t) exp i9G{x)) . (29) 

Since 

E, [exp {9G{x{aRit))))] > E, [exp (0G(x(aR(t))))l,^<t] 

= P^iaR < t} exp {9R) , 

we obtain the bound 

PxWr <t} < exp (7^TV(r)t) exp i9{E - R)) . 

As a consequence PxWr < t} ^ as R ^ oo and thus the Markov process 
x{t) is non-explosive. 

It follows that G{xR{t)) G{x{t)) almost surely as i? — > oo, so by the 
Fatou lemma we obtain from Eq. ( p9[ ) the bound Eq. (p6|). The bound Eq. 
( p7|) is obtained by noting that the left side is equal to 

PxWe{i+S) <t} < exp (70Tr(r)t) exp {-59E) , 



and this concludes the proof of Lemma 3.5. 

We have the following "tracking" estimates to the effect that the random 
path closely follows the deterministic one at least up to time tE for a set of 
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paths which have nearly fuh measure. We set Ax{t) = x{t,uj) — x^etit) = 
{Ar{t), Ap{t), Aq{t)) with both x{t) and x^etit) having initial condition x. 
Let 

S{x, E, t) = {x(-) ; G{x) = E and sup G{x{s)) < 2E} . 

0<s<t 



By Lemma gj, P{S{x, E, t)} > 1 - exp {-feTr{T)t - BE). 

Proposition 3.7 There exist constants Eq < oo and c > such that for 
paths x{t,uj) G S{x,E,tE) with tE = £'i/'=2/"i/2^ g^^^ E > Eq we have 



sup 

0<t<tE 



l|A<z( 
l|Ap(t)|| 
IIArf 



< c sup \\y/2^Tu 

0<t<tE 



E'-2 
_! 1 



V 



1 



(30) 



Proof: We write differential equations for Ax{t) again assuming both the 
random and deterministic paths start at the same point x with energy 
G{x) = E. These equations can be written in the somewhat symbolic form: 



dAq = Apdt , 

dAp = {0{E^-^/^^)Aq- K^Ar^dt, 
dAr = {-jAr + AAp) dt + y^^fTdu; 



(31) 



The 0{E^ 2/^2 -J coefficient refers to the difference between forces, —'VqV{-) 
evaluated at x{t) and x^etit)', we have that G{x{t)) < 2E, so that VqV{q) — 
VgV{qdet) = 0{d'^V)Aq = 0{E^-^/''^)Aq. For later purposes we pick a 
constant c' so large that 



p = p{x) = c E 



> supY^ sup 

i j {q:V{q)<2E} 



5V(g) 



dqidqj 



for all sufficiently large E. 

In order to estimate the solutions of Eqs. ([3T|), we consider the 3x3 
matrix which bounds the coefficients in this system, and which is given by 



M 




(32) 



We have the following estimate on powers of M; For AX(o) = (0,0,1)^, 
we set AX("^) = M™AXW. For a = max(l,7 + A), we obtain AX^^) < 
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a(0, 1, l)"^ AX(2) < a2(l, 1, 1)"^ and for m > 3, 




( 



< a Z 



m-2 

P 2 

m — 1 

p— 

171-2 



where the inequahties are componentwise. From this we obtain the bound 

(33) 




< 



i(ai)2ev^2ai 



l + at + i(at)2ev^2ai 



If < t < tE "we have < \/c'- Then the exponentials in the above 
equation are bounded, and 



(34) 



for some constant c. 

Returning now to the original differential equation system Eg. (|3l|) , we 
write this equation in the usual integral equation form: 





/ 

























v 






[ 






Ap(s) 

-v,y(g(s, cj)) ds + Vgy(gdet(s)) - r^Ar(s) 

-7Ar(s) + AAp(s) 



+ 



(35) 



From this we obtain the bound 



/ l|Ag(t)|| 
l|Ap(t)|| 
V l|Ar(t)|| 



< 



l|Ag(t)|| 
M| ||Ap(t)|| I ds + 
|Ar( 




where M is the matrix given by Eq.(|3^), and uj^ 
Note that the solution of the integral equation 



sup,<. ||V2^u;(i)ll- 



AX(t) = C dsMAX[s) + 



/ 







(36) 
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is AX{t) = exp(tM)(0,0,a;max)^- We can solve both Eq.(|5D and Eq.(|3|) 
by iteration. Let Axm{s), AXmis) denote the respective m*^ iterates (with 
Axo(s) = {0,0,^/2^uJ{s)f, and AXo{s) = (0, 0, w^nax)^, < s < t^j). The 
AXm's are monotone increasing in m. Then it is easy to see that 

/ „ 

^ ' < AX^{t) < AX{t), 



V 



\Aqm 

\Apm 

AV'fYi 



for each iterate. By Eqs. 
Eq. @ follows. I 



1), and the definition of p the conclusion 



As a consequence of Theorem 3.3 and Proposition 3.7 we obtain 



Corollary 3.8 Let Q.[E) = with a < l/k2 and assume that w{t) is such 
that supq<j<j^ ||v^27TLij(t)|| < Q,{E) and x{-,uj) £ S{x, E,tE)- Then there 
are constants c > and Eq < oo such that all paths x{t,w) with initial 
condition x with G{x) = E > Eq satisfy the bound 



{s)ds > cE" 



(37) 



Remark 3.9 For large energy E, paths not satisfying the hypotheses of the 
corollary have measure bounded by 



P,{ sup \\,/2^uj\\>n{E)} + I'{S{x,E,tEf} 

0<S<tE 




n{Ef 



+ exp {9{jTr{T)tE - E)) 



where a and b are constants which depend only on the dimension of to. Here 
we have used the reflection principle to estimate the first probability and 
Eq. ( p7D and the definition of S to estimate the second probability. For E 
large enough, the second term is small relative to the first. 



Proof. It is convenient to introduce the L -norm on functions on [0,t], \\j \\t = 
(/o t^s) • By Theorem U, there are constants Ei and ci such that 

for E > El the deterministic paths 2;det('S) satisfy the bound 



kdetllL = fj rl^t{s)ds>ciE^2 ^ 
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By Proposition 3.7, there are constants E2 and C2 such that ||Ar(s)|| < 
C2^{E), uniformly in s, < s < t^;, and uniformly in x with G{x) > E2. So 
we have 

Ikllt^ > llrdetlk^ - ||Ar||t^ > fciE^2 2j -c2n{E) f^^2 2 j . 

But the last term is 0{E'^''^/'^^^/'^^'^ ), which is of lower order than the first 
since a < l/k2, so the corollary follows, for an appropriate constant c and 
E sufficiently large. ■ 



3.3 Liapunov Function and Exponential Hitting Times 

With the estimates we prove now our main technical result. 

Theorem 3.10 Let s > and 9 < Oq = (max{ri, r„})"^ Then there 
are a compact set U = U{s, 6) and constants k = n{U, s,9) < 1 and L = 
L{U, s,6) < 00 such that 



T'exp{eG){x) < Kexp{0G){x) + Llu{x) 



(39) 



where Ijj is the indicator function of the set U. The constant k can be made 
arbitrarily small by choosing U large enough. 

Proof: For any compact set U and for any t, T^ exp (9G){x) is a bounded 
function, uniformly on [0,t]. So, in order to prove Eq.(^), we only have to 
prove that there exist a compact set U and k < 1 such that 

sup E^[exp(0(G(x(s))-G(x)))] < k < 1. 

Using Ito's Formula to compute G{x{s) — G{x) in terms of a stochastic 
integral we obtain 



[exp(0(G(x(s))-G(x)))] 



exp (07Tr(r)s)E^ 



exp ^—9 J dt + 9 J y^2jTrduj{t) 



.(40) 



For any 9 < 9q, we choose p > 1 such that 9p < 9q. Using Holder inequality 
we obtain, 



E, 



exp (^-9 jr^ ^^^^ Jq V^rduj{tyj 
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exp 



X exp 



exp —qO 



2 , QP^ 
7r at H 



X E, 



exp 



p 



.2n2 



[\^/2iTrf dt + Bp r 
Jo Jo 



E, 



exp 



r. 2 , QP^ 
dt 7r H 



y/2iTrduj{t) 



i/p 



dt{y/2YTr 



Here, in the next to last line, we have used the fact that the second factor 
is the expectation of a martingale (the integrand is non-anticipating) with 
expectation 1. Finally we obtain the bound 



< 



E, [exp(^(G(x(s)) 
exp (07Tr(r)s)E^ 



exp 



G(x)))] 
qO{l 



dt jr' 



1/9 



(41) 



In order to proceed we need to distinguish two cases according if 3//c2 — 
1/2 > or 3/k2 — 1/2 < (see Corollary 3^). In the first case we let Eq 



be defined by s = -E, 



l/fc2-l/2 




T. For E > Eq we break the expectation Eq. 



( ^l|) into two parts according to whether the paths satisfy the hypotheses 
of Corollary 



or not. For the first part we use Corollary 3.5 and that 
jQr^{s)ds >~/^^r2(s) > 0^3/^2-1/2. f^j. ^j^g second part we use estimate 
( |38| ) in Remark |3.9| on the probability of unlikely paths together with the 



fact that the exponential under the expectation in Eq. (41) is bounded by 
1. We obtain for all x with G{x) = E > Eq the bound 



E, [exp {e{G{x{s)) - G(x)))] < exp {ejTr(T)tE,) x 



exp 



[l-pOT„ 



)cE^2' 



+ a exp 



b'jtE 



i/g 



(42) 



Choosing the set U = {x; G{x) < Ei} with Ei large enough we can make 
the term in Eq. (^) as small as we want. 

If 3/A;2 — 1/2 < 0, for a given s and a given x with G{x) = E we 
split the time interval [0, s] into E^^^"^/^^ pieces [tj,tj+i], each one of size 
of order £'^/'^2-i/2g_ Yot the "good" paths, i.e., for the paths x{t) which 
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satisfy the hypotheses of Corollar y |3.8| on each time interval the 
tracking estimates of Proposition |3.7| imply that G{x{t)) = 0{E) for t in 
each interval. Applying Corollary 3.5 and using that G{x{tj)) = 0{E) we 
conclude that r2(s)ds is at least of order ^3/fc2-i/2 ^ ^i/2-i/fc2 = ^2/^2. 
The probability of the remaining paths can be estimated, using Eq. 
not to exceed 



1 — a exp 



1 _ 1 



The remainder of the argument is essentially as above, Eq. (42) and this 
concludes the proof of Theorem 3.1C. | 

The existence of the Liapunov function given by Eq.(|39[) can be inter- 
preted in terms of hitting times. Let tu be the time for the diffusion x{t) to 
hit the set U . 

Theorem 3.11 Assume that e < (max{Ti,r„})~^ For any (arbitrarily 
large) a > there exist a constant Eq = E{)[a) > such that for U = 
{x ; G{x) < Eq} and x G U'^ we have 



E,: \e 



< e" + e" 



l)exp ie{G{x) - Eo)) . 



(43) 



Proof: Let s = 1 and 9 < 9q he given, we set k = exp (— a)/2 and take U 
to be the set given by Thm. 3.1C| . Let Xn be the Markov chain defined by 
Xn = x{n) and Nu be the least integer such that Xj^^j E U . Then 



aN„ 



(44) 



so that to estimate the exponential hitting time, it suffices to estimate the 
exponential "step number". 

Using Chernov's inequality we obtain 



'PANu >n} = P,.{- ^(G(X,) - G(X,_i) < G{x) - Eq, Xj e U'} 



i=i 

n-l 



n 



,e(G{X,)-G{X,_i)) 
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E 



^0(G(X„)-G(X„_i) 



E, 



n-1 



< • • • < e 



e(G(x)-£;o) 



sup E,[e^(G(^^)-G(2^))] 



By Thm. 3.1C we have 



sup E,[e^(«(^i)-GW)] < K, 

and therefore we have geometric decay of P>„ = Fxi^jj > n} in n, P>„ < 
exp {9G{x) — Eq). Summing by parts we obtain 



E, 



:,aNu 



hm 

M~>oo 



J2 e'^^'PAru = n} 



n=l 



AI 



^a(n+l) _ an\ , a 



a(M+l) 



P~ 



>M 



.n=l 



which, together with Eq. gives Eq. (Esl). 



4 Accessibility and Strong Feller Property 

In this section we prove that the Markov process is strong Feller and more- 
over we show that it is strongly aperiodic in the sense that for all t > 0, all 
X & X and all open sets A C X we have Pt{x,A) > 0. Both results imply 
immediately that x{t) has at most one invariant measure: Since the process 
is strong Feller the invariant measure (if it exists) has a smooth density 
which is everywhere positive by the property of aperiodicity. Obviously no 
two different such measures can exist. 

The strong Feller property is an immediate consequence of the hypoel- 
liptic properties of the generator L of the diffusion. The result is an easy 
consequence of the estimates in [Q, H] , since there much stronger global hy- 
poelliptic estimates are proven (under stronger conditions on the potential 
We present here the argument for completeness. 
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The generator of the Markov process x{t) can be written in the form 



2d 



i=l 



If the Lie algebra generated by the set of commutators 



(45) 



has rank dim(X) at every point x £ X, then the Markov process has a C°° 
law. In particular it is strong Feller. This is a consequence of Hormander 
Theorem |11, |l^ or it can be proved directly using Malliavin Calculus de- 
veloped by Malliavin, Bismut, Stroock and others (see e.g. [|19|). 



Proposition 4.1 // H2 holds then the generator L given by Eq. ^T^ ) sat 

isfies the rank condition (45). 

Proof: This is a straightforward computation. The vector fields Xi, i = 
1, - ■ ■2d gives d (j) , i = l,n, j = 1, ■ ■ ■ , d. The commutators 



d 



1 



(j) 



Pi 



d 



u ) , Xq 
1 



1 -'^o 



yield the vector fields d o) and d (j) . Further 

Pi 91 



djj) , Xq 



, d (j)d (I) 

1=1 qY' q\' 



P2 



If is strictly convex, this yields d (j) while in the general case we need 
to consider further the commutators 



d 



9 u„ 
I1 



1=1 ^^(^'"'^gf" 



P2 



E 



Qm+ljj{2) 



,-1 d u,) ■ ■ ■ d u^)d (i) P2 

The condition H3 means that we can write d a) as a linear combination of 

P2 

these commutators for every x £ X. The other basis elements of the tangent 
space are obtained inductively following the same procedure. ■ 
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We now prove the strong aperiodicity of the process x{t). This is based 



on the support theorem of Stroock and Varadhan [27|. The support of the 
diffusion process x{t) with initial condition x on the time interval [0, t], is by 
definition the smallest closed subset Sx,t of C([0, i]) such that Px[x{t,uj) E 
Sx,t] = 1- The support can be studied using the associated control system, 
i.e., the ordinary differential equation where the white noise uj{t) is replaced 
by a control u{t) G L^([0,r]): For our problem we have the control system 

Q = P, 



P 



r = {—jr + Ap) + u , (46) 

and we denote x^it) the solution of this control system with initial con- 
dition X and control u. The support theorem asserts that the support of 
the diffusion Sx,t is the closure of the set {xu ', u G L^([0,t])}. As a conse- 
quence suppPf(x,-), the support of the transition probabilities is equal to 
the closure of the set of accessible points {y ; 3m S -L^([0,t]) s.t. Xu{t) = y}. 

Proposition 4.2 // condition HI holds then for all t > 0, all x ^ X 

suvvPtix,-) = X. (47) 

Proof: This result is proved in ||^] under the additional condition that the 
interaction potential C/^'^^ is strictly convex, in particular V?7^^^ is a dif- 
feomorphism. Our condition HI implies that S/U^"^^ is surjective. We can 
choose an inverse g : R"' R*^ which is locally bounded. From this point 
the proof proceeds exactly as in Theorem 3.2 of Q and we will not repeat 
it here. 



5 Proof of Theorem 



The proof of Theorem 2T is a consequence of the theory linking the ergodic 
properties of Markov process with existence of Liapunov functions, a theory 
which has been developed over the past twenty years. The proof of these 
ergodic properties relies on the intuition that the compact set U together 
with a Liapunov function plays much the same role as an atom in, say, a 
countable state space Markov chain. The technical device to implement this 
idea was invented in ||l|, ^ , and is called splitting It consists in constructing 
a new Markov chain with state space Xq U Xi, where Xi are two copies 
of the original state space X. The new chain possesses an atom and has 
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a projection being the original chain. The ergodic properties of a chain 
with an atom are then analyzed by means of renewal theory and a coupling 
argument applied to the return times to the atom. A complete account of 
this theory for a discrete time Markov process is developed in the book of 
Meyn and Tweedie [^], from where the result needed here is taken (Chapter 
15). 

For a given s > consider the discrete time Markov chain Xj = x{js) 
with transition probabilities P{x,dy) = Ps{x,dy) and semigroup = T^". 
By the results of Section ^, the Markov chain is strongly aperiodic, i.e., 
P{x, A) > Q for any open set A and for any x and it is strong Feller. The 



exponential bound on the hitting time given in Theorem |3.11| implies in 
particular that E2,[t[/] is finite for all x G X and thus we have an invariant 
measure (for hypoelliptic diffusions this is established in [|l4|). By aperi- 
odicity and the strong Feller property, this invariant measure is unique. 
The following Theorem is proved in 

Theorem 5.1 // the Markov chain {Xj} is strong Feller and strongly ape- 
riodic and if there is a function W > 1, a compact set U and k < 1 and 
L < oo such that 

PW{x) < kW{x)+L1u{x), (48) 
then there exist constants r > 1 and R < oo such that, for any x, 

Y^r^Pix,-)- ^i\\w < RW{x), 



where the weighted variation norm \\ ■ \\w is defined in Eq. (14)- 



By Theorem 3.10| the assumptions of Theorem |5.1| are satisfied with W 



exp{9G) and 9 < (maxjTi, T^})"^. For the semigroup T* we note that we 
have the apriori estimate T*exp(0G)(x) < exp(76'Tr(r)i) exp(0G)(a;), cf. 
Lemma |3.5| which shows that T* is a bounded operator on L^(X) defined 
in Eq. (|l5|). Setting t = ns + u with < n < s, and using the invariance of 
/i one obtains 

\\T' - fi\\e < ||T"^ - MT'We < Rr-' , (49) 
for some f > 1 and i? < oo or equivalently 

r*||Pt(x,-) -Mllexp(eG) < Rexp{eG{x)) . 

As a consequence, for any s > 0, has 1 as a simple eigenvalue and the 
rest of the spectrum is contained in a disk of radius p < 1. The exponential 
decay of correlations in the stationary states follows from this. 
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Corollary 5.2 There exist constants R < oo and r > 1 such that for all f , 
g with f'^, g"^ £ Lf{X), we have 



J fT'gdfi- J fdfi J gdfi 



<R\\frJw/'r-^ 



Proof: If G Lg°, we have < \\f'^\\y^ exp{9G{x)/2) and similarly for 

g. Further if Eq.(|9|) holds with W = exp (60) it also holds for exp (90/2) 
and thus for some Ri < oo and ri > 1 we have 



T'gix)- I gdfi 



<i?ir-||,^||f exp^'^^(^) 



Therefore we obtain 



J fT'gdfi- Jfdfijgdfi < J \f{x)\ T'g{x) - j g d^L 

J exp {eG)df?j Rirl 



dfi 



To conclude we need to show that / exp {9G)d^ < oo. This follows from 
Eq.(|48|) which we rewrite as 



eexp(0G(x)) < exp {9G{x)) - P exp {9G{x)) + Llu{x) , 
with e = 1 — K. From this we obtain 

. N . . N 

e- J2 exp {9GiXk)) < - exp + ^uiXk) • (50) 



k=l 



N 



N 



k=l 



By the Law of Large Numbers the r.h.s of Eq.(|50|) converges to Lfj,{U) which 
is finite, and thus / exp {9G) dfi is finite, too. ■ 



This concludes the proof of Theorem 2.1. 
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